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Abstract. Our main goal is the comparative study of singularities of solutions to 
the systems of first order quasilinear PDEs and their perturbations containing higher 
I— I derivatives. The study is focused on the subclass of Hamiltonian PDEs with one 

Oh spatial dimension. For the systems of order one or two we describe the local structure 

of singularities of a generic solution to the unperturbed system near the point of 
"gradient catastrophe" in terms of standard objects of the classical singularity theory; 
we argue that their perturbed companions must be given by certain special solutions 
of Painleve equations and their generalizations. 
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1. Introduction 

In [23! developed a perturbative approach to classification of systems of integrable 
Hamiltonian partial differential equations (PDEs). The systems under consideration 
are written in the form of long wave expansion (also called small dispersion expansion) 

(1.1) Ut = A{u) + 52(u; u^, u^^) + Bsi^u; u^, u^^, u^^^) + .... 
Here 

u = [u^{x,t), . . . ,u^{x,t)Y 

is the unknown vector function, the entries of the coefficient matrix A{vl) are smooth 
on some domain X> C M"; the characteristic roots Ai(u), . . . , A„(u) will be assumed to 
be pairiwise distinct 

(1.2) det (A(u) - A ■ 1) = 0, Ai(u)^Aj(u) for z ^ j, Vu G P. 

The terms i?2, -B3, • • • , -Bfc, ... of the expansioiij^ in the rhs of (1.1) are polynomials in 
the jet coordinates 

Ux = (m^;, ■ ■ ■ -iU^ , Mxx = [y-xx-i • • • -I '^xx) -I • • • 

graded homogeneous of the degrees 2, 3, fc, 
(1.3) 

deg — = m, m > 0, deg = 0, z = 1, . . . , n. 

ox''^ 

The coefficients of these polynomials are smooth functions on the same domain T). 

^In order to uniformize the notations the leading term could be redenoted as 

A(u) -Bi(u; u^;). 

However we will not follow this system of notations. We emphasize that the terms of degree zero 



(thus, depending only on u but not on the jets) must not appear in the rhs of ( 1.1 1 
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The system (1.1) can be considered as a perturbation of the first order quasihnear 
system 

(1.4) Ut = A{u)n.^ 

when considering slow varying solutions, i.e., working with vector functions u{x) chang- 
ing by 1 on the spatial scale of order L 3> 1. The magnitude |u(x)| will not be assumed 
to be small; the natural small parameter 

(1.5) i 

appears in the estimates for the derivatives 

(1.6) u^. ~ e, u^.^. ~ e^, . . . , u^^^ ~ e^ . . . . 
It is convenient to introduce slow variables by rescaling 

(1.7) X ^— ex, t ^ tt. 



The system (1.1) rewrites in the form 

(1.8) Ut = A{\i) Wx + t B2{u; u^, u^.^,.) + e^Bs{u; u^., u^;^;, u^xx) + .... 



In these notations the leading term (1.4) results from the limit e — * 0. 



For sufficiently small e one expects to see no major differences in the behaviour of 



solutions to the perturbed and unperturbed equations (1.8) and (1.4) within the regions 
where the x-derivatives are bounded. However the differences become quite serious near 
the critical point (also called the point of gradient catastrophe) where derivatives of the 
solution to the unperturbed equation tend to infinity. 

Although the case of small viscosity perturbations has been well studied and un- 
derstood (see [5] and references therein), the critical behaviour of solutions to general 



conservative perturbations (1.8) to our best knowledge has not been investigated in a 
systematic way before the author's paper |T7] (see the papers [291 Ell [361 EH EH] for 
the study of various particular cases). 

The main goal of this paper is to extend the technique of ^7] to Hamiltonian per- 
turbations of systems of PDEs. A serious difference from the scalar case treated in 
[T7] can be seen already in the case of second order systems considered in the present 
paper. For the scalar case any Hamiltonian perturbation preserves integrability up to 
a high order of the small parameter e. For systems this is not the case: for generic 
perturbation of a second order quasilinear system]^ integrability is destroyed already at 
the order e^. 

Due to this difference we begin our study of critical behaviour in perturbed sys- 
tems of Hamiltonian PDEs with considering integrable perturbations only. Our idea 
is that the local structure of the critical behaviour should depend only on the order 
e^ perturbation near a weak singularity of a hyperbolic system on the order e^ per- 
turbation near a weak singularity of an elliptic system. (We postpone to a separate 
publication the study of singularities on the boundary between domains of ellipticity 
and hyperbolicity.) Because of this it suffices to have integrability only at low orders 



^Recall that any second order quasilinear PDE (1.4 1 can be integrated by the method of 
characteristics. 
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of the perturbative expansion. We give a precise meaning to such an "approximate 
integrabihty" , explain the basic tools involved (extension of infinitesimal symmetries, 
D-operator, quasitriviality and string equation), and formulate the Universality Con- 
jectures about the behaviour of a generic solution to the perturbed Hamiltonian system 
near the point of gradient catastrophe of the unperturbed solution. We argue that, up 
to shifts, Galilean transformations and rescalings this behaviour essentially does not 
depend on the choice of solution (provided certain genericity assumptions hold true). 
Moreover, this behaviour near the critical point is given by the following formulae: 

• near a critical point within the domain of hyperbolicity the solution, after an affine 
transformation of the independent variables (x,t) (a;+,x_) and a nonlinear change 
of the dependent variables {u,v) ^ (r+,r„) the solution has the following asymptotic 
representation 

(1.9) 

where U = U{X; T) is a particular real smooth for a// X G M solution to the fourth 
order analogue of the Painleve-I equation 



1.10) X = TU- 



-U^ + — iU'^ + 2U U") + —U'^ 
6 24^ ' 240 



depending on the real parameter T (the so-called Pf equation). Here a, b, 
some constants that depend on the choice of the equation and the solution. For the 
original variables u, v the above asymptotic behaviour implies 

(1.11) 

with some constants A, /i. 

The solution U{X,T) was completely characterized by T.Claeys and M.Vanlessen 
in [9] via the isomonodromy deformations technique (see p^) where a long standing 
problem of existence of a pole-free real solution to Pf has been settlecj^ In the very 
recent paper [8] by T.Claeys and T.Grava the Universality Conjecture formulated in 
[17] was rigorously established, with the help of Deift - Zhou steepest descent method 
[13], for the case of critical behaviour of solutions to the KdV equation with analytic 
initial data. For the systems the somewhat weaker asymptotics ( |1.11[ ) was proved in 
[To] for a particular class of solution to the hierarchy of Toda equations with analytic 
initial data. These are the solutions appearing in the study of Hermitean matrix inte- 
grals. They depend on one arbitrary analytic function that determines the probability 



3n 



The equation (1.10 1 appeared in [B] (for the particular vakie of the parameter T = 0) in the study 



of the double scaling limit for the matrix model with the multicritical index m — 3. It was observed 



that the generic solution to (1.10) blows up at some point of real hne; the conjecture about existence 



of a unique smooth solution has been formulated. Our original conjecture formulated in [T7] says 



that the solution to (1.10 1 smooth on the entire real axis exists and is unique for any real T. The 



uniqueness part of this conjecture remains open. 
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distribution; the general solution to Toda hierarchy equations should depend on two 
arbitrary functions. The critical behaviour (1.11) appears, according to |T^, in the 



large asymptotic behaviour of orthogonal polynomials near the so-called singular 
edge points. The small dispersion parameter in this case is e = 1/A^. 

• near a critical point within the domain of ellipticity a suitable complex combination 
w = w{u, v) of the solution u{x, t), v{x, t) as function of an affine complex combination 
z of the independent variables behaves as follows 



w c^w° + ae^^^Wo (e 



-4/5 



+ O (e^/^) 



;i.i2) 



z = ax + bt + zo, a, 6 G C. 
Here Wo = Wo{Z) is the so-called tritronquee solution to the Painleve-I equation (P/) 

W" = 6W^-Z. 

This solution discovered by Boutroux in [4j is uniquely characterized by the following 
asymptotic behaviour 



W{Z) 



, , , , 47r 

-, \Z\ oo, I argZ| < — . 

In particular it has no poles for sufficiently large in the above sector. 

In [TS] it was conjectured that the tritronquee solution has no poles inside the sector 

4:71 



;i.i3) 



arg Z\ < 



Although a rigorous proof of this conjecture remains an open problem, numerical evi- 
dences found in [18j support the conjectur. 



Importance of absence of poles in the sector (1.13) has the following explanation. It 



turns out that the complex parameters a and 6 in (1.12) are such that for sufficiently 
small \t\ and e > the argument 

ax + bt + zq 



of the tritronquee solution belongs to the sector (1.13) for all real x. So the Universality 



Conjecture (1.12) makes sense only if no poles occur within the sector (1.13) 



The Universality Conjecture (1.12) appeared first in the recent paper [TH] for the case 



of focusing nonlinear Schrodinger (NLS) equation. In this paper we will explain that 
the same critical behaviour arises also for other Hamiltonian perturbations of second 
order quasilinear PDEs of elliptic type. 

The paper is organized as follows. In Section 2 we recall some basic facts about sec- 
ond order quasilinear PDEs, their Hamiltonian structures, solutions, conserved quan- 
tities and symmetries. In Section 3 we give a very simple local classification of weak 
singularities of solutions, for both hyperbolic and elliptic cases. An attempt to extend 
this local classification of singularities to quasilinear hyperbolic systems of order higher 
than 2 is also given; only the case of integrable hyperbolic PDEs (the so-called semi- 
hamiltonian systems in the terminology of [58]), or systemes rich according to |15] is 
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considered. In Section 4 the main tools for the study of Hamiltonian perturbations are 
developed. Applications of these techniques to the study of concrete examples is given 
in Section 5. 



Acknowledgments. This work is partially supported by European Science Foun- 
dation Programme "Methods of Integrable Systems, Geometry, Applied Mathematics" 
(MISGAM), Marie Curie RTN "European Network in Geometry, Mathematical Physics 
and Applications" (ENIGMA), and by Italian Ministry of Universities and Researches 
(MUR) research grant PRIN 2006 "Geometric methods in the theory of nonlinear waves 
and their applications" . 



2. Nonlinear wave equation: Hamiltonian structure, first integrals 

and solutions 



We now proceed to considering the two-component systems of Hamiltonian PDEs. 



In this section we start with the first order quasilinear systems (1.4). Moreover, we 
restrict ourselves with considering an important subclass of such systems associated 
with nonlinear wave equation 

(2.1) Uu-dlP{u) = Q 



for a given smooth function P{u). The equation (2.1) is linear for a quadratic function 
P{u)] we assume therefore that 

P"\u) ^ 0. 



2.1. Main examples of nonlinear wave equations. The equation (2.1) of course 
is of interest per se; it also arises in the study of dispersionless limits of various PDEs 
of higher order. In particular, for P{u) = — | one obtains the dispersionless limit 

(2.2) utt + {uu^)^ = Q 
of Boussinesq equation 

(2.3) Utt + {u u^)^ + u^.^.^.^ = 0. 



For P{u) = e" (2.1) reduces to the long- wave limit 

(2.4) Utt = 
of Toda equations 

(2.5) iln = e""+i^"" - e""~""-i 

More generally, one-dimensional system of particles with neighboring interaction 

(2.6) H = J2lpl + Piqn-qn-i) 

with the potential P{u) (generalized Fermi - Pasta - Ulam system) yields, after inter- 
polation 

quit) - qn-i(t) = w{en,et) 

(2.7) 

Pn(t) = v{en,et) 
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and substitution 
(2.8) 

the following system 

Ut = Vx 

(2.9) 



u 



1 - e-^^= 



w 



vt = e 



-1 



P' 



- 1 



^2 



-U 



P' 



1 - e 



-u 



d^P\u)^—[2P"{u)u^^^ + AP"'{u)u^u^^ + P'''{u)ul]+0{e' 



The above formulae are understood as formal power series in e: 



W 



''^ k>l ^ ' 



U 



(k) 



(2.10) 



u = w + 



A;>1 



Bk are the Bernoulli numbers. In the limit e — (2.9) reduces to the nonlinear wave 
equation written as a system 



Ut = fa 



(2.11) 



Vt = d^P'{u). 

Further examples will be considered in Section [5j 



2.2. Hamiltonian formulation of nonlinear wave equation. The system (2.11) 
can be written in the Hamiltonian form 

Ut = 



(2.12) 

with the Hamiltonian 
(2.13) 



vt = d. 



' 6v{x) 

6H 
Su{x) 



H 



-v^ + P(u) 



dx. 



The associated Poisson bracket reads 

(2.14) {u{x),v{y)] = 6'{x-y) 

The Poisson bracket of two local functionals 



dx, G 



dx 
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is again a local functional given by 

' 6F 



(2.15) 
where 



{F,G} 



^ 6G 6F ^ 6G 



6u{x) ^ 6v{x) 6v{x) ^ 6u{x) 



dx 



SF 
Su{x) 
5F 



SO 



Eu9 



6v{x) 



are the variational derivatives, 

E„ : = 



(2.16) 



d_ 

du 

F - ^ 

J^t) •— 7^ 
OV 



6u{x) 
ov[x) 



o^ - h 



d_ 



d'' 



d 



are the Euler - Lagrange operators. 

Recall that the densities of the local functionals are considered modulo total x- 
derivatives. So, the Poisson bracket (2.15) vanishes iff the integrand in the rhs is a 
total x-derivative. The following important criterion (see, e.g., [H]) gives a necessary 
and sufficient condition for a differential polynomial to be a total x-derivative of another 
differential polynomial. 

Theorem 2.1. Let 

M 

k,l>0 i,j 

be a differential polynomial with coefficients Aij{u,v) smooth on a disk D C Here 

i = {ii,...,ik), j = {ji,---,ji), l<ii,---,ik,ju---,ji<N 

are multiindices of the summation, M, N are two nonnegative integers. Then there 
exists another differential polynomial h of the same form such that 



iff 



5H 



m>0 



U 



(m+1) 



d 



+ v 



(m+l)_ 



d 



5u{x) 



6H 

6v{x) 



0, H 



hdx. 
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Corollary 2.2. Two local junctionals F = J f dx and G = J gdx commute with 
respect to the Poisson bracket (2.14) iff 



(2.17) 



E„ 



5F ^ 5G 5F ^ 5G 



6u{x) 6v{x) 6v{x) 6u{x) 



6F ^ 6G 

Or. 



6F ^ 6G 



6u{x) 6v{x) 6v{x) 6u{x) 







0. 



Alternatively (|2.1|) can be recast into the form of Euler - Lagrange equations 

'1 



5S = 0, S 




:0t - V{<P, 



dx dt 



with 



P(u) = V'(u) 



by introducing a potential 
The canonical Poisson bracket 



u. 



{(f){x),(f)t{y)} = S{x-y) 



yields d^lil). 



2.3. Domain of hyperbolicity and Riemann invariants. From (2.11) it readily 
follows that the nonlinear wave equation is hyperbolic on the domain of convexity of 

PM, 



(2.18) 



{u,v) G 



such that P"(u) > 



and elliptic when P{u) becomes concave. On the domain of hyperbolicity (2.11) can 
be reduced to the diagonal form 



(2.19) dtr± = ±y/P"{u)dxr± 
introducing Riemann invariants 

(2.20) r±=v±Q{u), where Q'{u) = ^P"{u). 

We also remind the following property of solutions written in the Riemann invariants: 
introducing the characteristic variables 

(2.21) x± = x± ^P"{u)t 
the function r^/r_ does not depend on 



(2.22) 



dx. 



Dxa 



0. 
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2.4. Solutions of nonlinear wave equations. Let us first describe the local struc- 
ture of general solution to the system (2.11). 

Lemma 2.3. Let f = f{u,v) be an arbitrary solution to the linear PDF 

(2.23) fun = P"{u)U. 

Let {uo,Vq) be a point such that 

(2.24) [fuv{uo,Vo)f - P"{uo) [U{uo,Vo)f ^ 0. 
Then the system of equations 

X = fu{u,v) 

(2.25) 

t = fviu,v) 

defines a unique pair of functions {u{x,t),v{x,t)) near the point 

ixo,to) := ifu{uo,Vo), fy{uo,Vo)) 

such that 

u{xo, to) = Mo, v{xo, to) = Vo. 



The functions {u{x,t),v{x,t)) solve ( |2.11 ). Their x- derivatives at the point {xo,to) do 
not vanish simultaneously. Conversely, any solution to (2.11) satisfying the condition 

(2-26) K-i'-M^a (.0.0)^0 

can be locally obtained by the above procedure starting from a suitable solution f{u,v) 
to the linear PDE (|2^. 

This is a reformulation of the classical method of characteristics (also called hodo- 
graph transformation). 



Remark 2.4. The Riemann invariants (2.20) are also the characteristic variables for 
the linear PDE (2.23), i.e., in these variables the PDE for the function f = f{r^,r^) 
reads 



(2.27) 



dr+dr^ 8 (P")3/2 ydr^ Or 



0. 



Remark 2.5. The exceptional solutions to (2.11) that cannot be written in the form 
(2.25) can be represented in the form 

V = Q{u) + c or V = —Q{u) + c 

for a constant c where u = u{x,t) satisfies the scalar quasilinear PDE 

ut = ^/P"{u)u^ or ut = -^/P"{u) u^. 

On these solutions one of the Riemann invariants takes constant values. 
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2.5. Conserved quantities and infinitesimal symmetries of nonlinear wave 



equation. An alternative interpretation of solutions to the linear PDE (2.23) is given 
by 

Lemma 2.6. Given a function f = f{u,v), consider the Hamiltonian 



(2.28) 



f{u, v) dx. 



This Hamiltonian commutes with the Hamiltonian (2.13) of nonlinear wave equation 



iff the function f satisfies (2.23) 



Proof. The Poisson bracket (2.14) of Hj and H reads 



{Hf,H} 



6u{x) 6v{x) 6v{x) 6u{x) 
The integrand is a total a;-derivative iff 

du{fu) = d,{f,P"{u)). 



dx = [fuVr, + f^P"{u)u^] dx. 



This gives (2.23). 



□ 



Recall that the differential polynomial / = f{u,v;Ux,Vx,...) is called conserved 
quantity (or density of a first integral) of the PDE ( |2.11[ ) if the following identity holds 



(2.29) 



dtf = d^g, dt = Y^ 



m>0 



for some differential polynomial g. It is understood that in (2.29) one has to substitute 



from (2.11) Ut ^ fx, Vt i— > P"{u)u^ 



Corollary 2.7. The function f = f{u,v) is a conserved quantity for the nonlinear 
wave equation (2.11) iff it satisfies the linear PDE (2.23) . 



Proof. For a solution / to (2.23) one has 

dtf = d^g 

where g = g{u,v) is determined by a quadrature 

dg = P"{u)fvdu + fudv. 
The converse statement is straightforward. 

Corollary 2.8. The Hamiltonian PDE 

Us = d^f^{u,v) 

(2.30) 

Vs = d^fu{u,v) 

commutes with the wave equation 



□ 



{ut). 



[U 



sJt ' 



(vt), 



^s)t 



if the function f{u,v) satisfies the linear PDE (2.23) 
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Proof. Indeed, the flow (2.30) is generated by the Hamiltonian Hj: 

6Hf 



a 



6v{x) 
6u{x) 



As it is well known from Hamiltonian mechanics, commutativity of Hamiltonians im- 
plies commutativity of the flows. □ 

2.6. Commutative Lie algebra of infinitesimal symmetries of nonlinear wave 



equation. We conclude that the solutions to the linear PDE (2.23) are associated with 



infinitesimal symmetries of the nonlinear wave equation (2.1). A somewhat stronger 
statement says that all Hamiltonians of the form Hf commute pairwise, i.e., the Lie 
algebra of symmetries is commutative. 



Lemma 2.9. Given two arbitrary solutions f, g to the linear PDE (2.23), the Hamil- 
tonians 



H 



f 



commute: 



f{u,v)dx and Hg 



{Hf,Hg} = 0. 



g{u, v) dx 



We leave the proof of this lemma as an exercise for the reader. 

For a given P{u) one obtains therefore an infinite family of commuting Hamiltoni- 



ans parametrized by solutions to the linear PDE (2.23). The associated Hamiltonian 
systems commute pairwise. 



Remark 2.10. All equations (2.30) of the commutative family have the same Riemann 



invariants (2.20).- 
(2.31) 



dsr± = {huv ± \/P"{u)K^) dxr±. 
Example 2.11. Given a real number k 7^ 0, ±1, choose 



(2.32) 

The function 
(2.33) 



Pu) 



2 



u 



k{k + 1 



satisfies the linear equation (2.23). The associated Hamiltonian flow reads 

ut = d^{uv) 

(2.34) 
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These equations coincide (after the change of sign of the time variable) with equations 
of motion of one- dimensional isentropic gas with the equation of state 



Here u = p (the mass density), v the velocity. Therefore the equations (2.34) of gas 
dynamic can be considered as an infinitesimal symmetry of the nonlinear wave equation 



u 



K-l 



,K-2„,2 



Utt = -U^x + U U^. 

K — i 

Example 2.12. For the exceptional value k = 1 take 

(2.35) P{u) = Tu{\ogu~l). 
The choice 

(2.36) f = l{u^^uv^) 



of a solution to the linear PDE (2.23) produces the dispersionless limit 

(2.37) Ut + {uv)x = Q 

of the focusing /defocusing nonlinear Schrddinger equation 

(2.38) iiJt + \^.x±\i^?i^ = ^ 
written in the coordinates 



I.e., 

ut + {uv)x = 

Vt + VVxT Ux = 



u 2 u 



So, the dispersionless NLS is an infinitesimal symmetry of the following nonlinear wave 
equation 

Remark 2.13. In [H] a hierarchy of the first order commuting quasilinear PDEs has 
been associated with an arbitrary solution F = F{v^, . . . ,v"') of equations of associa- 
tivity 

d^F d^F d^F d^F 

nPl — riP'^ j A h. 1 — ^ „ 



dv^dv^dvP dv'idv'^dv'- dv'-dvWvP dvidv^'dv 

(2.39) 

d^F 
dv^dv^dv'j 
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Here Vj^^ are entries of a constant symmetric nondegenerate matrix, rjij are entries of 
the inverse matrix. The quasihomogeneity axiom has not been used for the construction 
of the hierarchy. For n = 2 the solution F = F{v, u) must have the following form 



for some function For 



^'(u) = P(u) 



this family of commuting PDEs coincides with the one associated with solutions to 



(2.23). The Hamiltonians of the hierarchy o/ [H] are solutions to (2.23) polynomial in 



V. 



Before proceeding to the discussion of the critical phenomena let us extend the 



method of characteristics described in Lemma |2.3| to solving any equation of the com- 
mutative family 



(2.40) 



with the Hamiltonian 



Vg Oxhr^ 



-f^/i = y h{u,v)dx, huu = P"{u)h^.o- 



Lemma 2.14. Let f = f{u,v) be any solution to ( |2.23 ) satisfying the condition 

- fLhlY - P"{uo) - Chi)' + 

for some (uo,fo)- Then the system 

X + shuv{u,v) = fu{u,v) 



(2.41) 



sh^^{u,v) = fv{u,v) 



u{xo,So)=Uo, v{xq,So)=Vo 

defines a unique pair of functions {u{x, s),v{x, s)) near (xq, sq) determined from 

Xo + So huv{Uo, Vo) = fu{Uo, Vq) 



These functions solve the system (2.40). Conversely, any solution to (2.40) satisfying 



certain nondegeneracy conditions can be obtained by this procedure. 
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3. Singularities of solutions to nonlinear wave equation 



In this section we develop a local classification of singularities of solutions to the 
nonlinear wave equation. Actually, the technique based on the implicit function theo- 



rem works equally well also for any PDE of the form (2.40) from the commuting family 



of infinitesimal symmetries of the wave equation. So, our study of singularities will be 



developed for the solutions to the general system (2.40). 



The functions (m(x, s),v{x, s)) determined by the system (2.41 ) are smooth provided 
the conditions of the implicit function theorem hold true. At the points (xc, Sc) where 



the implicit function theorem fails to be applicable the solution to (2.40) has a weak 
singularity, i.e., there exist the limits 



lim 



M X, S) 



lim f (x, s) 



but the derivatives Ux, Vx blow up. The main goal of this Section is to describe the 
local behaviour of a solution near the generic weak singularity. The precise definition 
of genericity will be given below. 

Let us first subdivide weak singularities into three classes. 

Type I. The limiting value {uc, Vc) belongs to the domain of hyperbolicity of the 



system (2.40), i.e. 
(3.1) 



P"{ 



Ur 



> 0. 



Type II. The limiting value {uc, Vc) belongs to the domain of ellipticity of the system 



(2.40), i.e. 



(3.2) 



P"(ur) < 0. 



Type III. The limiting value {uc,Vc) belongs to the boundary between the domains 
of hyperbolicity and ellipticity. 



(3.3) 



P"{ 



0. 



Observe that the linear PDE (2.23) belongs to the hyperbolic/elliptic type near the 



critical point {uc,Vc) of Type I/II, while the type of this PDE changes near a generic 
(i.e., P"'{uc) 7^ 0) critical point of Type III. 

In this Section we will study the local behaviour of the solution near the weak 
singularity following the above classification. Only Type I and Type II singularities 
will be considered; the study of Type III singularities is postponed to a subsequent 
publication. 

Let us do first some general simplification of notations. First, since the partial 



derivative satisfies the linear PDE (2.23) if the function h does so, it suffices to 



study singularities of implicit functions ^(x, s) and f (x, s) determined by the system 



(3.4) 



X + shu{u,v) = fu{u,v) 



shy{u,v) = fv{u,v) 
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where both functions h and / satisfy (2.23). Next, the shifts 

X X — Xc, S ^— > S — Sc 

can be absorbed by a redefinition 

f ^ f + XcU + Sjl 

of the function /. So, without loss of generahty we may assume that Xc = Sc = 0. 

We will also adopt the following system of notations: the values of the functions 
f{u,v), h{u,v) and of their derivatives fu{u,v), fy{u,v) etc. at the critical point 
(tic,fc) will be denoted /°, /°, /° etc. We will also denote 

Po:=P{u,), P',:=P'{u,) 

etc. 

According to the above conventions we will have 

(3.5) /° = /° = 0. 

Violation of the implicit function theorem conditions in our notations reads 

(3.6) - det ( : ' \- ^™ : ' J;™ ) = (/d ^ - p^ {cf = o. 

3.1. Critical points of Type I and Whitney singularities. As it follows from the 



representation (2.25) of solutions, the inverse map to 

(x,t) 1-^ {u{x,t),v{x,t)) 

is well defined also on a neighborhood of the critical point. Moreover, it is given by 
the gradient 

(3.7) {u,v) ^ {f^{u,v)J^{u,v)) = (x,t) 



of a function satisfying the linear PDE (2.23) that, in the case of Type I singularity 



at (uc, Vc) belongs to the hyperbolic type. The points of blow-up of the derivatives of 
the solution {u{x,t),v{x,t)) are exactly the points where the Jacobian of the gradient 
map ( |3.7 ) vanishes. 



Classification of singularities of smooth maps of plane to plane obtained in the cele- 
brated paper of H.Whitney |49| was one of the first successes of the singularity theory. 
Recall that Whitney obtained the complete list of normal forms of smooth maps 

(M,0) ^ (M,0) 
{u,v) ^ {x,t) 

stable with respect to small perturbations. He proved that, up to diffeomorphisms 

(x, t) ^— (£, t) 
(m, v) ^— s> (m, v) 

these maps belong to one of the following three normal forms: 
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Type Wl 
(3.8) 

Type W2 
(3.9) 

Type W3 
(3.10) 



X = V 



t = u 



X = v'^ 



t = u 



X = uv 

6 



t = u 



For the Wl case the inverse functions u = u{x,t), v = v{x,t) are well defined and 
smooth. We will show below that also in the general situation the solutions to the 



shape (3.9). 



quasilinear system (2.40) near a generic Type I singularity locally has the universal 



Another natural question is about the structure of solutions near the first point of 
singularity. We say that (xc,tc) is the first singularity of the solution to the nonlinear 
wave equation if the functions u{x, t) and v{x, t) are smooth for sufficiently small |x— xd 
and for t < tc- We will see that the local behaviour of a generic solution near a first 
point singularity, within the hyperbolicity domain, is described by the W3 singularity 
on the Whitney list. This sounds like being almost obvious; an additional point is that 
only affine transformation of independent variables x, t will be allowed near the critical 
point. 



Let proceed now to the study of singularities of solutions to the general system (2.40 ) 



Due to the assumption Pq > the jacobian (3.6) factorizes 



(f y 

\J uv / 



yjuv^y^ojvvj \Juv y ^0 Jvv j 



The first of the genericity assumption says that, at the critical point only one of the 
factors vanishes. Without loss of generality we may assume that 



(3.11) 





uv 



Ji 



II /O 

vv 



0. 



fuv + \/P^fvv 7^ 



Definition 3.1. We say that the critical point (|3.11|) of Type 1 is generic if: 
(3.12) 
and 

(3-13) — \/ Pq f^y^ + /^t, ^ p„ 7^ 



pill 

^0 
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or 



hut 



(3.14) 



J uvv V J vvv ' Ji 



TDlll 







uv 



pill 

"P^hl + K^J^^O and 



J uvvv 



pill 

J WW ' J vt 



wv ^ pii 



+ 





w 



32 (p^')5/2 



The formula 
(3.15) 



dxr± 



fuv ^ fw ^ {^^v ^ ^/^f^vvj 

for the x-derivatives of Riemann invariants clarifies the meaning of the first genericity 



assumption (3.12): at the critical point (xc, Sc) = (0,0) only one Riemann invariant 



blows up (under our assumptions r_ does so) while another one remains smooth. 



We now show that the case (3.13) is incompatible with the assumption that (xc, Sc) 



(0, 0) is locally the first critical point. 



Lemma 3.2. Let {xc, Sc) = (0,0) be a critical point of Type I satisfying (3.11), (3.12) 



such that the functions u{x,s), v{x,s) determined by (3.4) along with normalization 
u{xc,Sc) = Uc, v{xc,Sc) = Vc are smooth for s < Sc for sufficiently small \x — xd, 
\s — Sr\- Then 



(3.16) 



f 

J uv 



fpTi fO I ^0 
Q J vvv ' J v\ 



TDlll 
w ^pu 



0. 



The meaning of the constraint (|3.16|) is the following: the graph of the Riemann 



invariant r_(x,Sc) must have an inflection point at x 
Proof. For simplicity of calculations we will assume that 



Under this assumption one can solve the system (3.4) in the form 

fuhy fvhu 



s{u, v) 



k 

hy'' 



X = x{u, v) 



Expanding these functions in Taylor series in r+, r_ at u 
following expressions for the linear combination^ 

(3.17) x± = x+fe± 



Vr we obtain the 



Clearly these correspond to the characteristics of the hyperbolic PDE. 
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of independent variables 



"-v \ Juvv ' V ^0 Jvvv i ^Jvv \ "■uv ' V -^0 "'vv 



fl + 0{\w\^) 



(3.18) 



-pill 
8 J"' 



4 1/: 



fO / p/// 







/p7i-0 _|_ ^0 -^0 



f2 



0{\w\^) 



where we introduced the notations 



(3.19) 



r± = r± — r_|_ 



for the shifted Riemann invariants. Assuming (3.13) do a rescahng 

x± ^ k x±, f+ I— >■ /cf+, f„ (—> k^^'^f_ 
with a small positive /c. After simple calculations one obtains 



(3.20) 



where 



1/2N 



2;+ 



2a;_ 



J 1 







/T?" f _|_ f -^0 

^ ./ vvv ' J vv 4 p" 



(observe that 7^ due to the genericity assumption (3.12)). We obtain Whitney 
normal form of type 14^2. Therefore the solution to the implicit function system (3.4) 
near the critical point is defined on a half plane with the boundary a;_ = passing 
through the critical point. Clearly such a behaviour is incompatible with the assump- 
tions of the lemma. □ 



Theorem 3.3. Under the genericity assumption (3.14) the local structure of the solu 



tion to the implicit function equations is described by the Whitney normal form of type 
W3. 



Proof We have to include the degree 3 terms in the Taylor expansions of x{u, v) 
s{u,v). The rescaling 



(3.21) 



x+ 



kx-, 



r_ I — > 
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yields 
(3.22) 

X. = -%Aof+f^ + iBofI + Oik'/'] 
h^, 6 



where Aq and Bq are defined in (3.14). The theorem is proved. □ 



3.2. Critical points of Type II and elliptic umbilic catastrophe. Recall |17] 
that elliptic umbilic catastrophe can be represented as the singularity a.t X = Y = 
of the function 

(3.23) F{X,Y) = ^X' - XY^. 

It can also be identified (see [2]) as one of the real forms of the codimension 2 singularity 
of type D4. The following interpretation will be useful in sequel. Let $(X, F) be a 
harmonic function of two variables, 

"^xx + = 0. 
Consider singularities of the gradient map 

(3.24) (x,y) ^ ($x(x,y),<i>y(x,r)) = iu,v) 

i.e. the points where the Jacobian of ( |3.24[ ) vanishes, 

/ dU/dX dUdY \ 

det = - [<i>j,^ + <^j,y] = 0. 

\ dV/dX dV/dY J 

As it follows from the above formula for the Jacobian, at this point all the second 
derivatives of $ vanish. So, the 3-jet of the function must be a cubic harmonic poly- 



nomial. It must coincide with (3.23) up to a linear change of variables 



The gradient map of (3.23) is given by the following system 

U =X^-Y^ 

(3.25) 

V =-2XY 
In the complex coordinates 

W = U-iV, Z = X + iY 

this gives the quadratic map 

(3.26) W = 

We will now show that the Type II singularity of a generic solution to the quasilinear 



system (2.40) is locally described by the elliptic umbilic catastrophe. 

As above we assume that the singular point is xc = sc = 0; it will also be assumed, 
as above, that 
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Denote u = u — uc, v = v — vc- Consider a solution to (2.40) determined by the system 



(2.41) near the critical point. 



Lemma 3.4. At the critical point all second derivatives of f vanish: 
(3-27) = = = 0. 



Proof follows from (3.6) due to negativity of Pq. 



Definition 3.5. The Type II critical points is called generic if the 3-jet of f{u,v) does 
not vanish at this point. 



Lemma 3.6. After the rescaling 



X 

s 
u 

V 



k X 
k s 
k^/^u 
k'/^v 



at the limit k 0+ the solution to the system (3.6) near a generic Type II critical 



point {x = Xc = 0, s = Sc = 0,u = uc, v = vc) tends to the solution of the equation 



(3.28) 



z := x + shl + icosh'^^, w:=v + icoU, oq := /l^ + i Cq 



where 



Co 



Observe that, due to (2.23) and (3.27) at the critical point the third derivatives of 



the function f{u, v) satisfy 

J uuu oJuvV J uuv •! vvV 

So, at the generic Type II critical point we have necessarily ao 7^ 0. 



Proof of Lemma is similar to that of Lemma 3.2 Expanding the left hand sides of 



eqs. (3.4) in Taylor series and using equation (2.23) valid for both functions / and h 



one obtains 



x + shl + s {-clhlu + hl,v) + - [^(-c^M^ + v^) _ 2clhl,,u v] 



sh', + s (hlu + hlv) + - [hl,i-cy + v') + 2hl^u v] 
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After the rescaling and division by k one arrives at the system 

x + shl =\ [fl^{-clu^ + v^) - 2clfl,uv] + O (A:V2) 

sK =\ [/L(-cgn2 + v') + 2fl,u v] + {k'/') . 
This converges to the quadratic equation (3.28). The Lemma is proved. 

We see that, after a hnear change of independent variables the local structure of 
the solution near a generic Type II singularity is described by the bifurcation diagram 
(3.26) of the elliptic umbilic catastrophe. 

3.3. Generalization to the multicomponent case. Let us consider an n-component 
system of PDEs written in the diagonal form 

(3.29) dtUi = ai{u)dxUi, i = 1,. . .,n, 

n > 3. Thus the dependent variables ui, . . . , m„ are the Riemann invariants for the 
system (3.29); the characteristic velocities ai(u), a„(u) will be assumed to be 
smooth functions taking pairwise distinct values, 

aj(u) 7^ aj(u), i j for all u eV. 

Recall the following definition 



Definition 3.7. The system (3.29) is caWed semihamiltonian if the coefficients satisfy 
the following differential equations 

(3.30) dk ( r^^l = ^ 



ai — Qj J \ai — Ofc 

for all triples of pairwise distinct indices i, j, k. 

Here 

. ._ A .-rl . 

we will also use similar notations for higher derivatives, i.e. 

Oijfc '■= djdkai, Ai jki := djdkdiAi 

etc. 

The meaning of this definition is clarified by the following theorem 

Theorem 3.8. (i) Let Ai(u.), . . . , An{u) be an arbitrary solution to the overdetermined 
system of PDEs 

(3.31) A,, = {A, - A,) , 1 + 3. 
Then the system of PDEs 

(3.32) dsUi = Ai{\i)d^,Ui, z = l,...,n 



commutes with (3.29), 

(3.33) dsdtUi = dtdsUi 
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(ii) For a semihamiltonian system (3.29) the general solution to (3.32) depends on n 
arbitrary functions of one variable. 



(iii) Given a solution to the system (3.31) such that 
Ai(uo) = ■■ 



A„(uo) = 0, JJ/li,,(uo) ^0 



for some Uq G T), then the the system of equations 

(3.34) X = ttiivi) t + Aiiyi) , i = l,...,n 

defines a unique smooth vector function u{x,t) for sufficiently small \x\, \t\ such that 
u(0,0) = uo. This function satisfies the system [3.29] . Conversely, any solution u{x,t) 



to the semihamiltonian system (3.29) satisfying 



9,ui(0,0) ^0,...,9,M„(0,0) ^0 
can be locally obtained by the above procedure. 



Remark 3.9. For a semihamiltonian system (3.29) there also exists an infinite family 



of conserved quantities h{u) such that the following identity holds true on the solutions 



to (3.29) 



(3.35) dMu) = dj{n) 

for some function /(u). The conserved quantities are determined from the following 
overdetermined system of PDEs 



(3.36) 



didjh 



dih + 



djh, 



The solutions to this system also depend on n arbitrary functions of one variable [ibid]. 
The conserved quantities h{u.) were called entropies in [45j; the systems (3.29) possess- 
ing a rich families of entropies (i.e., families depending on n arbitrary functions of 
one variable) were called systemes rich in [45j. This class coincides with the class of 
semihamiltonian systems. 

Let us describe the local structure of singularities of solutions to the semihamiltonian 
systems near the point of gradient catastrophe. Without loss of generality we may 
assume that the gradient catastrophe takes place at the point Xq = 0, to = 0. 

At the point of catastrophe the assumptions of the implicit function theorem applied 



to the system of equations (3.34) fail to hold true. Explicitly this means that 

n 



(3.37) 



det {aij{u)t + Aij{u)) 



at the point of catastrophe. 

Definition 3.10. The point x = 
called generic if at this point 

(3.38) 



x=xo, t=to, u=uo 







i=l 



A, 



Xq = 0, t = to = 0, u = Uq of gradient catastrophe is 
= 0, A,-,(uo)^0 for jj^i, 
(uo) = 0, Ai^iii{uo) ^ 0. 
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The first part of tliis definition says that only z-th Riemann invariant breaks down at 
the point of catastrophe. The assumption Aj_jj(uo) = is not a restriction; it says that 
the graph of the i-th Riemann invariant Ui{x, 0) has an inflection with vertical tangent 
at the point of catastrophe. As we will see the assumption Ai^m^uo) ^ guarantees 
that the inflection point does not degenerate. 

Let us assume that a; = xq = 0, t = to = 0, u = uq = is a generic point of gradient 
catastrophe with i = n. We will denote 



a° := ai(uo), := ^^^^(uo 



._ 



etc. Observe that 



Alj = for all i ^ j 



due to the system (3.31) 



Perform a linear change in the space of independent variables putting 



(3.39) 

Theorem 3.11. After the reseating 



Z'i • — 3^ Ct^t ^ Z — 1 ^ . . . ^ Tl^ . 



(3.40) 



t ^ 

Un I ^ k ^ Un. 



the system (3.31) converges to 



n 



(3.41) 



Proof. Spelling the equations (3.30) out one obtains 



/„ Ar,\ — aj — ttk {ai — ajYai^kO'kj — {a^i — o^kYo^i^jaj^k 

\oAZ) aijk — 77 — (ii,jO'i,k + - 



'''^ {tti - aj){ai - ttk) 



{ai - aj){ai - ak){aj - au) 



for all triples j, of pairwise distinct indices. Using these equations along with the 
system (3.31) yields 
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Lemma 3.12. At the point of catastrophe one has 



0, -4^,„i = for m ^ n, i ^ m, n 



(3.43) 



0, A 





m,nni 



0, a: 





m,nij 



n aO _ n 



for m ^ n, i,j, k ^ m,n, k distinct 



^Li = 0, ^n,i,- = 0, for i,j j^n, i ^ j 



a: 



- =0 A° 

n,nm '-n,nij ^' n. 



0, A" 



ijfe 



for i,j,k^n, i,j,k distinct 



Corollary 3.13. The Taylor expansions of the equations (3.34) near the point of catas- 
trophe read 



9 2^ ^0 _' + Mm 2^ 



/ a. 



m,i ,0 



4^ II- 



/,0 _ ^0 ^'■m,m"'i 



1 a*^ 

I _ /lO 2 _|_ "m,?! .0 



rn,mtt i 



(3.44) + ] A^yy^-fniUi + ^m,mmm'^m ~^ ^ (l^l ) 



4^ 7/ 

■^m,m""m 



X — a„t = t 



i 



(3.45) 



It is understood that the summation indices i, j etc. in the sums denoted are all 
distinct from n and (in (3.44)) from m. 



Applying the rescaling (3.40) to the equations (3.44), (3.45) one completes the proof 
of the Theorem. □ 
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4. Perturbations of nonlinear wave equations 

4.1. Deformation problem-I: main definitions and the programme of re- 
search. Given a system of the first order quasihnear PDEs 

(4.1) Ut = A(u)u,, u= (ui(x,t),...,u"(x,t)) 
admitting a Hamiltonian description 

(4.2) ut = {u{x),Ho}o, Ho = jho{u)dx 

with respect to a Poisson bracket of hydrodynamic type written in the flat coordinates 
in the form 

(4.3) {u\x),u^{y)}o = rf^'{x - y), r]^' = r]'^ = const, det {r]'^) ^ 0, 



we say that a system of the form (1.8) is a Hamiltonian deformation of (4.1) if it can 
be represented in the form 

(4.4) ut = {u{x),H} 

with a perturbed Hamiltonian 

H = Ho + eHi + e^H2 + ... 

(4.5) 

Hk= hk{u;u^,...,u'^''^)dx, k>l 



deg/ifc(u; u^, . . . jU^'')) = k 
and a perturbed Poisson bracket 

{u^{x),u^{y)} = {u\x), u\y)}o + e {u\x), u^{y)}, + e^u\x), u^{y)h + • • • 

(4.6) 

fc+i 

{u%x),u^{y)}k = J2 (u(^); u.(a;), . . . , u^^\x)) 6^'~^+'\x -y), k>l 

s=0 

deg4^'Ju;u,,...,u(^)) =s. 

In the above formulae the notations like /(u; u^, . . . , u*^^)) are used for differential 
polynomials 

/(u; u,, . . . , u(^')) G C°°(M)[u,, . . . , u G M, 

i.e., for polynomial functions on the jet bundle J^{M). The degree of a differential 
polynomial is defined by the rule 

deg< = l, deg<^ = 2, i = l,...,n. 



The sums in (4.5) and (4.6) are allowed to be infinite. However every term of the e- 



expansions of (4.4) is a well-defined differential polynomial. Antisymmetry and Jacobi 



identity for the Poisson bracket are understood as identities for formal power series 
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in e. The delta-function symbol can be spelled out in the following way. Introduce a 
matrix of linear differential operators depending on the parameter e 



IT^ = IV^ + elV^ + e^n^ + . . . 



(4.7) 



k+l 



5^Ay,(u;u.,...,uW)e^+\ k>l. 



s=0 



One clearly has 
(4.8) 



{u'{x),u^{y)} = WUix-y). 



Then the perturbed Hamiltonian system reads 

6H 



(4.9) 



* 6uHx) 



m>0 k+l=m 



The assumptions about the structure of the perturbative expansions in (4.5) and (4.6) 
imply that the terms of the expansion ( |1.8[ ) 

6Hi 



(u; u., . . . u(-+i)) = 5^ n 



m > 0, 



fc+i=m 

are differential polynomials of the degree m + 1. 

The class of Hamiltonian deformations is invariant with respect to (special) Miura- 
type transformations of the dependent variables 

(4.10) u^^u = u + ^e'=Ffc(u;u^,...,u('=)), degF^ (u; u^, . . . , u^'^)) = /c. 

k>l 

The transformation of the Hamiltonian is defined by the direct substitution; the Poisson 
bracket is transformed by the rule 



(4.11) 



{u\x),u^{y)} = W^Six-y) 



where L and L"^ are respectively the operator of linearization of the transformation 



(|4.10|) and its formal adjoint: 

J-ii 



(4.12) 



The transformations of the form (4.10 ) form a group. Two Hamiltonian deformations 
of the quasilinear system (4.1) related by a transformation (4.10) are called equivalent. 
In particular the Hamiltonian deformation is called trivial if it is equivalent to the 
unperturbed system (4.1). 
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Example 4.1. The focusing/ dejocusing nonlinear Schrddinger equation 

(4.13) zeV'* + 2e'V^..±l^r^ = 

written in the coordinates 

I ,|2 e f tpa: i'x 



2i ytjj ijj 



takes the form of a Hamiltonian perturbation of the system (2.37): 
(4.14) ut + {uv), = 

^ ,2 



4 \ u 2 
The Poisson bracket gets no perturbation: 

{u{x),v{y)} = 6'{x - y), {u{x),u{y)} = {v{x),v{y)} = 0. 
The perturbed Hamiltonian reads 



(4.15) H = Ho + e^H2, Ho = J ^{u^ T uv^) dx, = - J 



— dx. 
8u 



Example 4.2. The generalized FPU system with the Hamiltonian (2.6) in the trans- 
lation invariant variables 

(4.16) Vn=Pn 

and interpolation 

Wn = w{en,et) 
Vn = v{en,et) 

yields the following equations for the functions w = w{x,t), v = v{x,t) 

1 11 

wt{x) = - [v{x) - v{x - e)] = Va:- -ev^^ + -(?v^^^ + . . . 
e 2 6 



(4.17) 



v,[x) = - [P'{w{x + e)) - P'iw{x))] = d,P'{w) + le dlP'{w) + dlP'{w) + .... 
e 2 6 



This system is a Hamiltonian perturbation of the nonlinear wave equation (2.11). The 
Poisson brackets of the variables w{x) , v{x) have a nontrivial dependence on e. Indeed, 
the brackets of the discrete variables read 

other brackets vanish. After interpolation and division by e one obtains the perturbed 
bracket 



{w{x),v{y)} = ^ [6{x - y) - 5{x ~ y - e)] 



(4.18) 



5\x-y)-\e5"{x-y) + \eH"'{x-y) + 
2 6 
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One can reduce the Poisson bracket (4.18) to the standard form 



{u{x),v{y)} = 6'{x - y) 



by the Miura-type transformation (2.8). Indeed, introducing the shift operator 

A = e'^^ 



rewrite the bracket ( 4.18[ ) as follows 



{w{x),v{y)} = - (1 - A 1) 6{x-y). 



The substitution (2.8) represented in the form 



u 



w 



1 - A-i 
implies 

{u{x),v{y)} = ^ {w{x),v{y)} = dj{x-y). 

In the new variables the deformed system is described by the perturbed Hamiltonian: 

ut = 



vt = d. 



6v{x) 

5H 

5u{x) 



H 



hdx 



v^{x) + P {w{x) — w{x — e)) 



dx 



(4.19) 



h=^v^ + P{u) 



f2 f4 

P"(u)ul + 



[8P"{u)ul-P'''{u)ut]+0{e' 



24 ' ' ^ 5760 

(modulo inessential total derivatives). The series in e in this case are all infinite. 
In fact the following general statement about deformations of the Poisson bracket 



(4.3) holds true. 



Theorem 4.3. 1). Any deformation (4.5) of the Poisson bracket (4.3) is trivial. 



2). Any special Miura-type transformation (4.10) canonical with respect to the bracket 
{,} = {, }o 

{u\x),u^{y)} = {u\x),u^{y)} 
is the time-e shift generated by a Hamiltonian 

^2 



u^=u' + e {u'{x), $} + - {K(x), $}, $} + 



(4.20) 



if dx 



99 = ^eVfc(u;u^,...,u('')), degv9fc(u;u^,...,u('')) = A;, A; > 

A;>0 
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with some differential polynomials ipo{u), (pi{u; u^), . . . . 



Proof follows from triviality in positive degrees of Poisson cohomology of the bracket 



(4.3) established in [25], [12] (see details in 



Thanks to the Theorem one reduces the study of Hamiltonian deformations of the 
Hamiltonian quasilinear equations to the study of systems written in the standard form 



(4.21) 



6H 



just perturbing the Hamiltonian, 
(4.22) 

iyW = y"/iW(u;u„...,u('=))(ix, deg/iW(u;u,,...,u('=)) = A;, A; > 

with some differential polynomials h^^^{u), h^^^{u; u^.), /^[^^(u; u^, u^x), ■ ■ ■ ■ Here 

is the Hamiltonian of the unperturbed equation. The equivalence of perturbed Hamil- 



tonians is established by the canonical transformations (4.20) 



(4.23) 



H^H = H + e{H,<^} + - {{H, $}, $} 



In the subsequent sections we will discuss the main problems of the deformation 
theory of Hamiltonian PDEs, namely. 

Problem 1: describe first integrals/infinitesimal symmetries of a given deformed 
Hamiltonian system by studying the procedure of extension of integrals/symmetries of 
the unperturbed system to integrals/symmetries of the perturbed system. 

Problem 2: classify integrable Hamiltonian perturbations for which all first in- 
tegrals/infinitesimal symmetries of the unperturbed system can be extended to first 
integrals/infinitesimal symmetries of the perturbed one. 

Problem 3: comparative study of the properties of solutions of the perturbed versus 
unperturbed systems. In a more specific way, characterize the perturbations for which 
the solutions to the perturbed Cauchy problem, in a suitable class of initial data, 
exists for the times beyond the moment of gradient catastrophe for the solution to 
the unperturbed Cauchy problem. For such a class of perturbed Hamiltonian systems 
establish the universal properties of solutions near the point of gradient catastrophe of 
the unperturbed system, that is, those properties independent from the choice of the 
initial data. A more ambitious project is in establishing universality independently on 
the choice of Hamiltonian perturbation. 
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4.2. Deformation problem-II: extension of infinitesimal symmetries and D- 
operator. Let 

ffprt = [°l + e ffW + + . . . 

,4.24) //-'l^//,>'>(u;u.....u.^.)... degft>^-(u; , . . . u<«) ^ .>0 

be a perturbation of the Hamiltonian of the nonhnear wave equation. We know that, 



given a solution / = f(u, v) to the linear PDE (2.23) the Hamiltonians iJ' ^ and 
(4.25) 
commute: 



H 



f — f f dx 



{i/[°l,iff} = 0. 
The goal is to construct a deformation 

(4.26) = j hf{u-\i^,...,M^^^)dx, deg/ifl(u;u^,...,u(*^)) = A;, fc>0 



for an arbitrary solution / of (2.23) such that 



{-ffpert, Hf} = 0. 

Definition 4.4. We say that the linear differential operator 



(4.27) 



DM=id, Z}W = y6fl^ (u;u,,...,u( . ^ . 
degfoi;! ... (u;u.,...,u('=)) = A;, A; > 1 



zs a D-operator for the perturbation (4.24) if for any solution f to (2.23) the Hamil- 
tonian 



(4.28) 

commutes with Hp^rt- 



Hf= Df dx 



In the formula (4.27) m{k) is some positive integer depending on k. It is easy to see 
that 



(4.29) 



m{k) 



'3k' 
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The summation is taken over all indices ii, im{k) from 1 to 2. As usual the 
coefficients bf^ ^ ^ (u; u^, . . . , u^'^)) are differential polynomials. As the D-operator 



makes sense only acting on the kernel of 

(4.30) dl - P"{u)dl 

the coefficients are not determined uniquely. One can choose a suitable normal form 
for this operator, e.g., keeping only derivatives of the form 

9™ and 

We will show that the D-operator is determined uniquely just from the commuta- 
tivity with iJpert, as it follows from the following 

Theorem 4.5. The D-operator, if exists, is determined uniquely up to a total x- 
derivative, i.e., if Di and D2 are two D-operators then 

D2f -DJ elmd^ \/feKei{dl-P"{u)dl). 

Proof follows from the following result saying that the nonlinear wave equation has no 
ffist integrals depending nontrivially on the derivatives. 

Lemma 4.6. Let P{u) be a potential such that P"{u) 7^ 0. Let 

^ = ^(u;u^.,...,u('=)) 

be a differential polynomial of the degree k > such that the functional 

G = j g{vL] u^.,. . .,VL^^^)dx 

commutes with H^^\ 

(4.31) {G,i7[°]} = 
Then g is a total x-derivative: 

g{u; u,, . . . , u^'^)) = d.J{n; u,, . . . , u^^"^)). 

The proof of this lemma is straightforward and will appear elsewhere. 
Corollary 4.7. Let the perturbed Hamiltonian (4.24) admit a D-operator. Then for 



any pair of solutions f , g to (2.23) the Hamiltonians Hf and Hg commute: 
{Hf, Hg} = 0, Hf = j Df dx, Hg = j Dg dx 

(4.32) 

fuu P i^ujfvvy guu P i^u^gw 

Proof From Jacobi identity it follows that 

{{Hf,Hg},Hpert} = 0. 

Assume the Poisson bracket {Hf,Hg} to have the form 
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for some k > 1. The functional -F^*^' must commute with i^'^l. Due to the Lemma it 
must vanish. □ 



Remark 4.8. One can consider a more general class of the first order quasilinear 
Hamiltonian systems written in the form 

ut = dxK 

(4.33) vt = d^hu 

for a given smooth function h = h{u,v). The conservation laws f = f{u,v) for this 
system are determined from the linear PDF 

(4.34) h^yfun huufyy. 

The associated Hamiltonians 



H 



f 



fdx 



commute pairwise. We say that the system (4.33) is regulaiH if it has no nontrivial 



local conservation laws depending on higher derivatives (cf. Lemma 4-6) ■ It would be 
interesting to characterize all regular quasilinear Hamiltonian systems. 

We end this Section with considering a simple example of Hamiltonian perturbations 
of the scalar transport equation 

Ut + UUx = 

classified in [T7j. The Hamiltonian structure of the transport equation reads 



Ut + {u{x),H^''^}=Ut + d, 
H^^i = - I u^dx 







6u{x) 

{u{x),u{y)} = 6'{x - y). 



All perturbations of the Hamiltonian H^^^ up to the order are parameterized by two 
arbitrary functions of one variable c = c{u) and p = p{u) as follows: 



pert 



U 

~6 



24 



c{u) ul + e^p{u) u. 



dx. 



This perturbation turns out to be integrable with the order approximation 

{Hf,Hg} = O (e^) , Hf = j D,Jdx, Hg = j D,,,gdx, V/ = f{u), g = g{u). 
The associated D-operator has the form 



D,J = f--cf"ul + e 



(4.35) 



24 



+ 



;2 y(4) 

480 



+ 



/(6) p'fW p;(5) 



+ 



+ 



1152 1152 3456 6 6 
^Similarly to the definition of regular elements in Lie algebras. 



ut 



o («•) , 
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the function / in this case is arbitrary. 



4.3. Deformation problem-Ill: an integrability test. Existence of a D-operator 
for a given perturbed Hamiltonian gives rise to a simple integrability test, as we will now 
explain for the example of the generalized FPU system. The perturbed Hamiltonian 
reads 



(4.36) 



pert 



Piu) 



24 



P"{u) ul 



dx + O (e^) 



(see (4.19) above). Under what condition for the potential P{u) a D-operator ex- 
ists within the approximation? We will show that this is possible only under the 
constraint 

(4.37) P"P^^ = {P 



iii\2 



To derive (4.37) let us look for the perturbed Hamiltonian density with indetermined 
coefficients 

Hf = j hfdx, {Hf,H^,,,} = 0{e^) 

hf = f + ehf{u; u^) + e'^hf{u; u^.,u^^). 

The first correction must be linear in Ux, Vx- Adding if necessary a total x-derivative 
one can reduce to the study of the first perturbation of the form 

= p{u, v) Vx 

for some function p = p{u, v). Computation of the bracket gives 

{Hf, ffpert} = e j Pu[vl- P"{u)ul] dx + O (e^) . 

It is easy to see that the integrand is never a total x-derivative unless Pu = 0, i.e. 
p = p{v). Hence hf must be itself a total x-derivative. 

Let us now consider the order 2 terms. Up to a total x-derivative they can be written 
in the form 



1 



(a(M, v)ul + 2b{u, v)uxVx + c{u, v)vl) . 



Computation of the order terms in the Poisson bracket gives the following functional: 

{Hf, if pert} = 



f S — f p"u + 

I I 12 "'XXX ~ 



f,P"'-bP"] Ux 



(cP"ux + bvx) Vx 



+ ^ {fvP'"" + l2P"a, - 24P"6,) ul-^ (a„ + 2P"c„) ulvx - \ (2a, + P" c,) Uxvl 



+ ^(c„-2y t;^! dx\0[e^) 



Denote / the integrand. One must verify the equations 

E„J = 0, E,/ = 0. 
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Vanishing of the coefficient of Uxxx in E^,/ gives 

C - —fv^ P". 

Next, vanishing of the coefficient of VxxVx in E^J imphes 

h = c 

Hence it exists a potential A = X{u,v) such that 

b = Au, c = A^,. 

After the substitution we obtain, collecting the coefficients of UxxUx, UxxVx, VxxUx, 
VxxVx, the following system of equations 

6P'"A, + P7™. = 

P"'Xu — P" \uu + P" ^vv = 

fi X — f,p" \ _i_ p" f — n 



It follows that 



A -J^f 
6 P'" 

A -JElf 

6 P'"''™"' 



Equating the mixed derivatives 
implies 



///^2 



fuvv • 



This proves (4.37). 



The equation (4.37) for the potential can be easily solved: dividing by P"'P" yields 

p'" pIV 



pn pin 

Integration gives 

P'" = cP". 

Hence 

P{u) = ke'^' + au + b 

for some constants a, b, c, k. 

We have proved that the only generalized FPU system passing the integrability test 
must have exponential potential. This gives Toda lattice. 

Remark 4.9. It would be interesting to compare the above arguments with the theory 
of small amplitude normal forms for finite FPU chains developed recently by A.Henrici 
and T.Kappeler [301 l3T]. 
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4.4. Deformation problem-IV: on classification of integrable perturbations. 

Let / = f{u,v) and g = g{u,v) be two first integrals of the nonlinear wave equation 
(2.11), i.e., two solutions to the linear PDE (2.23). The Hamiltonians 

Hf^ = J fdx, Hf = j gdx 

commute: 

We obtain a commutative Lie algebra of Hamiltonians parameterized by solutions to 
(2.23). The goal is to classify deformations of this algebra of the form (4.4) preserving 
the commutativity. 

So far the classification is available only within the approximation. Let us describe 
this preliminary result and outline the classification procedure. 



Theorem 4.10. Any 2-integrable perturbation of the nonlinear wave equation (2.11) is 
parameterized by two arbitrary functions of one variable p± = p±{r±). The associated 
D -operator reads 



Dp,rf = f + 



P i^Pufvvv ~l~ Pvfuvv^ ~^ 2^ Pijf'" 



( ff P'" 

I P Pvfvvv ~l~ Pufuvv ~l~ ^pii P^f^ 
+ {Pufvvv + Pvfuvv) vl] + C(e^), 



where 



P = p(m, v) = p+(r+) + p- (r_), 
the Riemann invariants r± = r±{u,v) are defined in (2.20). 



Proof Let / = f{u,v), g = g{u,v) be two arbitrary solutions to the linear PDE (2.23). 
Any (9(e)-perturbation of the commuting Hamiltonians and uf'^ can be written 
in the form 

Hf = Hf^ + e j p{u, v)vx dx + O (e^) , Hg = Hf + e j q{u, v)v^ dx + O (e^) 

for some functions p and q. From the commutativity 

{Hf,Hg} = 0{e') 
it follows, calculating the coefficient of u^x in Ji^i^f^ ^g}^ ^^^^ 

(4.38) ^ ^ 

Jvv Qvv 

Consider now a canonical transformation 

u ^ u + e{u(x), fsT} + O (e^) 
generated by the Hamiltonian 



K 



k{u, v) dx. 
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The Hamiltonians Hf and Hg will transform according to the following rule: 
Hf^Hf + e{Hf,K} + 0{e'), Hg ^ Hg + e{H^^\ K} + O {e') . 
Choosing the function k = k{u,v) satisfying 

kuu P {u)kyy — — — — 

Jvv 9vv 

one eliminates the terms linear in e from both perturbed Hamiltonians Hf and Hg. 

At the next step one has to consider the perturbation of order that can be chosen 
in the following form: 

Hf^Hf + '^j [aful + 2hfU,v, + Cfvl) dx + O (e^) 

Hg^H^g^ + '-^J {agul + 2bgu,v, + Cgvl) dx + O {e^) 

with some functions Oj = ai{u,v), bi = bi{u,v), Ci = Ci{u,v), i = f, g. From the 
commutativity 

{Hf,Hg} = 0{e') 
it follows, calculating the coefficient of v^xx in J^i-^f^ -^g}^ ^^^^ 

(a/ - CfP"{u)) g^y = {ag - CgP"{u)) f^^. 

Hence 

a/ = CfP"{u) + k{u,v)fyy 
ag = CgP"{u) + k{u, v)gyy 

for some function k{u,v). 

At the next step we substitute 

bf — 1"lfvvv ~l~ 1"2fuvv ~l~ l^sfvv ~l~ 1"4:fuv) Cf — Sifyyy -\- S2fuvv ~l~ ^sfvv ~l~ ^aJuv 
bg — figvvv ~l~ f2guvv ~l~ fzQvv ~l~ fiQuv-i Cg — Signyy + S2guvv ~l~ ^39vv ~l~ ^4guv 

with some coefficients n = ri{u,v), Sj = Si{u,v), i = 1, ... ,4 (the coefficients of the 
D-operator). The calculations similar to above yield the following expressions for the 
coefficients: 

p/// 

ri = P"py, r2 = Pu: rs = — p„, r4 = 

Sl = Pu, S2 = Pv, S4 = 

k = ^P"'p,-2s3P" 

where the function S3 = Ss{u,v) is arbitrary and the function p — p{u,v) satisfies the 
linear PDE 

pill 

Puu P Pvv 2P" 
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It is readily seen that this equation is equivalent to the D'Alembert equation in Riemann 
invariants: 

P'" d'^o 
-P"o -—0 =-4P"—-J— 
Pun ^ Pvv 2p// dr+dr_ ■ 

Finaly, we obtain the e^-correction of the Hamiltonian densities in the following form 



h 



f 



f + 



P i^Pufvvv ~l~ Pvfuvv) ~^ 2^ Pvfv 



"l~2 I P Pvfvvv ~l~ Pufu 



P' 



4P 



// Pu fv 



+ {Pufvvv + Pvfuvv) vl + S3 [vl - P"ul) fvv]+0 (e^) , 

a similar expression holds for hg replacing f ^ g. It is easy to see that the term 
containing the arbitrary function S3 can be eliminated by a canonical transformation 
generated by the Hamiltonian 

e 



with 



(u, v) determined from 



(j){u, v)vx dx 



S3- 



□ 



4.5. Quasitriviality and solutions to the perturbed nonlinear wave equations. 

Recall (see Section 4.1 above) that the e-perturbation of a first order quasilinear system 
is called trivial if it can be eliminated by a (special) Miura-type transformation (4.10). 
All terms in the e-expansion of this transformation are differential polynomials. 

Definition 4.11. We say that the perturbation is quasitrivial if it can be eliminated 
by a transformation of the form 

(4.39) u ^ li = u + 5^ (u; u,, . . . , u^^^)) , deg (u; u,, . . . , u("^W)) = k 

k>l 

where the terms for k > 1 are rational functions in the jet variables u^, u^x, 
. . . u^^^'^'^^ for some integer m{k). 

As it was proved in [19], all bihamiltonian perturbations are quasitrivial. Moreover, 
the order m{k) of the highest derivative in Fk can be estimated as 

'3k' 

(4.40) m{k) = — 

It would be interesting to extend the Quasitriviality Theorem of [19] to all Hamiltonian 
perturbation^ 

Existence of a quasitriviality transformation gives a universal formula for the per- 
turbative solutions to the deformed systems with very general classes of initial data. 



®The quasitriviality theorem was proved by S.-Q.Liu and Y.Zhang in [¥1] for a more general class 
of perturbations of scalar PDEs. In the recent paper |40| the quasitriviality technique was applied to 
the analysis of existence of a local Hamiltonian structure for scalar evolutionary PDEs. 
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Namely, substituting a solution uq = uo(s, t) of the unperturbed system into (4.39) one 
obtains a formal series in e satisfying the deformed PDE in all orders. In the examples 
of Section |5] below all the rational functions in the quasitriviality transformation have 
the denominators that are powers of 

(4.41) A = vl- P"{u)ul. 

In the Hamiltonian case the quasitriviality transformation can be represented as a 
time-e canonical transformation (see [U]). That means that there exists a Hamiltonian 

m>l 

(4.42) 

such that the canonical transformation 

(4.43) u ^ li = u + e{u(x), K} + -{{u{x),K}, K} + . . . 
transforms all unperturbed commuting Hamiltonians into the perturbed ones, 

(4.44) Hf = Jf{u)dx for any / satisfying fuu = P"iu)fvv 
The generating Hamiltonian K has to be found from the equation 

(4.45) Hf + e{Hf\K} + 'l{{Hf\K}, K} + ■ ■ ■ = Hj ^ Hf + e Hf + e'Hf + .... 

The Hamiltonians constructed in this way will automatically commute; the highly 
nontrivial set of cancellations ensures that the densities of these Hamiltonians up to 
total x-derivatives are 'polynomials in the jet coordinates 

The above perturbative expansions can be used for asymptotic integration of the 
deformed Hamiltonian systems for sufficiently small times (although rigorous analytic 
results justifying validity of these asymptotic expansions are still to be established). 
They clearly lose their validity when approaching the time of gradient catastrophe 
since the denominators become all equal to zero. In order to study the behaviour 
of the perturbed solutions near the critical point we will now describe an alternative 
procedure of constructing solutions to the perturbed Hamiltonian systems. Let us first 
prove the following 

Lemma 4.12. Let the perturbed Hamiltonian of the nonlinear wave equation have the 
form 



(4.46) i/pert= I \v^ + P{u) + Y,^^h^'\^-u,^^^,...M'\v^^^) 

L k>l 



dx 



(i.e., no explicit dependence on v in the higher corrections). The the D -operator, if 
exists, has v -independent coefficients. 
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Proof Indeed, under the above assumptions the coefficients of the commutativity equa- 
tion {Hf, i/pert} = for Hf = J Df dx 



E„ 




ou{x 
5K 







P"iu)u, + ed, 



pert 



6u{x) 



do not depend on v (see (2.17[)). Here 



k>l 



.u^'^Kv^'^^dx. 



Moreover, the shift 



V t-^ V + c 



maps to itself the space of solutions to the linear PDE (2.23). Therefore the same 
shift applied to the D-operator produces another D-operator. Due to the uniqueness 
theorem 4.5 the shifted operator coincides with D after adding, if necessary, a total 
x-derivative. □ 

Remark 4.13. The class of Hamiltonan perturbations considered in the Lemma cor- 
responds to the following subclass of perturbed nonlinear wave equations: 



(4.47) 



Utt 



P"{u)uxx + td^aiu] u^, ut, Uxx, Uxt, 



Given a solution h = h{u, v) to the PDE (2.23), consider a Hamiltonian perturbation 
of the quasilinear system 

ut = dxK 
vt = d^hu 

(we have changed notation s t for the time variable). Solutions to this system are 
parameterized in the form (2.40) by functions / = f{u,v) satisfying the same PDE 
(2.23) according to the method of characteristics. Let us assume that the perturbed 
Hamiltonian equation admits a f-independent D-operator. We will show that the 
solutions to the perturbed equation close to (2.40) can be obtained from a system of 
"infinite order ODEs" written in a variational form. This system of ODEs will be 
called string equation. 

Theorem 4.14. The solutions to the string equation 



x + t 



5Hy 



6H 



f 



6u{x) 5u{x) 



(4.48) 



6H^ 

6v(x) 



6H^ 

6v{x) 
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also solve the Hamiltonian equations 



(4.49) 



ut = 



vt = d. 



5Hh_ 
5v{x) 



where f = f{u,v) is another solution to fuu = P"{u)fvv, 



h'-= — 
dv ' 



Proof Since the coefficients of D-operator are f -independent, one has an identity 



(4.50) 

Due to obvious commutativity 



d 

—Dh = Dh'. 
ov 



ov 



^ F 

ov 



= 



m 



the assumptions of the Theorem from (4.50) one derives 

(4 51) d_m^^m^ d 6Hh ^ 6Hh> 

dv Su{x) 6u{x)' dv6v{x) 6v{x) 
Let us use the above identities in order to prove that the flow 



u. 



(4.52) 



td^ 



1 + td, 



6Hh' 
6v{x) 

6u{x) 



IS an 



infinitesimal symmetry of the system (4.49). Indeed, 



6Hf, 



{Ut)s = t {da:^-^, Hh'} + d, 

dv[x) 



d 5Hh 

" dv 6v{x] 



= t{{uix),Hh},Hh'} + d, 



6v{x) 



ivt)s = t{{vix),Hh},Hh,} + d,^^^^ 
Computation of the derivatives in the opposite order yields 

dv[x 



6H I 

{vs)t = dx^-/\ + {{vix), Hh'}, Hh}. 
du[x) 

The Jacobi identity and commutativity of the Hamiltonians {Hh,Hh'} = imply 

{{u{x),Hh},Hh'} = {{u{x),Hh'},Hh} 
{{vix),Hh},Hh'} = {{vix),Hh'},Hh}. 



42 



B. DUBROVIN 



So 



{U 



s)ti 



Js t- 



The commutativity remain valid adding to the symmetry (4.52) an arbitrary flow of 



the hierarchy, e.g., replacing (4.52) with 



Mo 



5Hy 

5v{x] 



-a 



5Hj_ 
5v{x) 



(4.53) 



5Hy 

5uix) 



6u{x) 



for some solution / to the linear PDE ( 2.23[ ). It is wellknown (see, e.g., [HI [37]) that 
the set of stationary points of an infinitesmal symmetry is invariant for the flow (4.49). 



Integrating one obtains the string equation (4.48). In the dispersion limit the string 



equation coincides with the equation (|2.41|) of the method of characteristics: 

fu 



X I h/ry^' 

th, 



4.6. Special Painleve-type transcendents and Universality Conjectures. We 

have seen in Section [3] above that locally near the point of catastrophe all generic 
solutions to a nonlinear wave equation and, more generally, to any Hamiltonian PDE 
commuting with the wave equation, have a very standard behaviour. Namely, if the 
catastrophe takes place in the domain of hyperbolicity then the local behaviour is 
described by the Whitney singularity equations. For the catastrophe in the domain of 
ellipticity the local behaviour is given in terms of elliptic umbilic singularity. In both 
cases the shape of the solution is locally universal, i.e., independent from the choice of 
the solution, up to a simple affine transformation of the independent variables and a 
nonlinear transformation of the dependent variables. 

The Universality Conjecture for the perturbed systems says that, similarly, the lead- 
ing term in the asymptotic expansion of the solution at e ^ near a critical point 
does not depend on the choice of a generic solution, up to simple transformations of 
the independent/dependent variables. The universal shapes that occur in such a local 
description are given in terms of certain very particular solutions to Painleve equations 
and their generalizations. Let us first describe these solutions. The description will be 
given separately for catastrophes taking place in the domain of hyperbolicity or domain 
of ellipticity. 

• Type I singularity. Consider the following 4th order ODE for the function U = 
U{X) depending on the parameter T 



(4.54) 



X = TU 



1 



6 



U^ + — (f/' + 2U U") + —U 



IV 



24 



240 



U' 



dU 
dX 



etc. 



The equation (1.10) can be considered as a higher order analogue of the Painleve-I 
equation. It is often identified as the Pf equation. As it was proved by T.Claeys and 
M.Vanlessen [9], this equation has a particular solution without singularities on the 
real axis X G M for all real values of the parameter T. Denote this solution U{X,T). 
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We refer to the paper [9] for the precise characterization of the particular solution in 
terms of the Stokes multipliers of the associated with (1.10) linear differential operator 
with polynomial coefficients. For large |X| the solution grows like 

U{X,T) ~ {-QXy/\ 

It is monotone decreasing for T << 0; for positive T it develops a zone of oscillations 
(see Fig.l). 



Fig.l: The solution to (1.10) for T > 



We are ready to formulate the Universality Conjecture describing the local behaviour 



of solutions to the perturbed Hamiltonian system (4.49) near the Type I critical point 



{xc,tc,Uc,Vc). Denote r± the Riemann invariants (2.20) of the unperturbed system. 
Let 

X± = {X - Xc) + C±{t - Q, C± = K{Uc,Vc) ± ^/P"{Uc)K{Uc,Vc) 

be the characteristic variables at the point of catastrophe. Let us assume that at 
the point of catastrophe of the unperturbed system the second Riemann invariant 
r_ breaks down, while the invariant r_ remains smooth (see (3.22) above). Then one 



should expect the following asymptotic behaviour of solutions to (4.49) near the critical 
point: 



6/7\ 



(4.55) 



r_(uc,yc) + a_e^/'^f/ (b^ 



X- 



for some constants a±, b±, c+ depending on the choice of solution and on the perturbed 



Hamiltonian. Here U = U{X, T) is the particular solution to the 4th order ODE (1.10) 
described above. 

In the next Section we will explain how to compute the constants for the particular 
examples of integrable systems. 
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• Type II singularity. Here we need to specify a particular solution W = Wq{Z) to 
the classical Painleve-I (P/) equation 

(4.56) Wzz = QW^-Z. 

The particular solution in question was discovered by Boutroux in the beginning of the 
last century It is characterized by the asymptotic behaviour 

/Z , , 47r , , 

argZ < — , Z — >• oo 
V 6 5 

(the principal branch of the square root is chosen). In particular the solution Wq{Z) 
is free of poles within the sector |argZ| < ^ for sufficiently large \Z\. Boutroux 
called it the tritronquee solution. A.Kapaev [M] gave a complete characterization of 
the tritronquee solution in terms of particular values of the Stokes multipliers of the 



linear differential operator with polynomial coefficients associated with (4.56). 




Fig.2: The real (left) and imaginary parts of the tritronquee solution to Painleve-I 

equation in the sector | axgZ\ < ^. 

In [T8] it was conjectured that all poles of the tritronquee solution to Pj are confined 
within the sector | arg Z\ > The Fig. 2 borrowed from shows the graphs of the 
real and imaginary parts of the tritronquee solution suggesting that, indeed there are 
no poles in the sector | argZ| < 

The Universality Conjecture for the Type II critical behaviour suggests the follow- 
ing structure of the leading term of the asymptotic behaviour of (complex) Riemann 
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invariants 

(4.57) r+ = r+(u,,v,) + P e'/' Wo (^) + O {e'/') , 

the dependence of the complex conjugate Riemann invariant r_ on the variable x_ is 
given by complex conjugation. Here Wo{Z) is the tritronquee solution to Pj, 

x+ = {x- X,) + (/i° + z v^^/^S) {t - tc), 

a, P are some complex numbers. These numbers are such that for any suffuiciently 
small \t\ and e > the line 

on the complex plane belongs entirely to the sector 

I arg Z\ < — 
5 

where the tritronquee solution Wq{Z) conjecturally has no poles. 

In the next Section we will consider in more details the most widely known examples 
of Boussinesq, Toda, NLS and Ablowitz - Ladik equations where are approach can be 
applied. For the case of Boussinesq equation we will give more detailed explanations 
that motivate our Universality Conjectures. We will see that, although these equations 
look so different, the local structure of critical behaviour in all these cases is given 
by the same formulae. The rigorous mathematical justification of such Universality is 
still to have been done; a more challenging problem is to establish Universality for a 
more general class of Hamiltonian perturbations, without restricting to only integrable 
perturbations. 



5. Main examples 



5.1. Boussinesq equation. As it was already explained in Section 2A_, the Boussi- 
nesq equation can be considered as a Hamiltonian perturbation of the nonlinear wave 
equation with the potential 



P{u) = --u\ 



The linear PDE (2.23) for this potential coincides with the Tricomi equation 



(5.1) 
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The domain of hyperbolicity is the half-plane m < 0. The D-operator transforming 
solutions of this equation to first integrals for Boussinesq hierarchy reads 

fvv ~l~ ^ fuvv ) l^x fvvv'^x'^x fuvv'^x 



(5.2) 



'6 



g ^ /4j; fuvv j ^xx g fu3v ~l~ '^fvvv) '^xx'^xx ^ ^ f'^^v'^xx 



1 2 2 1 

hvUxxVl - fu3vVxxUl - (31 Ai; " 20li^/6t, + 92u fu4v) U 

2 . 1 

-— m(13/5^ + 5ufu5v) ulVa: + — (H /n4t; " lOwVe?;) 

io iU 

+73 (7 /5„ + lOu fu5v) u^vl + fe^v^ 

15 D 



2 2 



Here and below we adopt the following notations for high order partial derivatives 



f4v 



9V 



u3v ■ 



du dv^ 



etc. In particular, the choice 
gives the Hamiltonian density of Boussinesq equation 



X 1 2 1 3 
•' 2 6 



1 1 
DbouJ = 2 ~ 6 ^ 2"^^' 

Let us proceed to the study of critical behaviour of solutions to an equation of the 
Boussinesq hierarchy 

5H 



(5.3) 

where 
(5.4) 



ut = 



vt = 



dv{x) 

SH 
5u{x) 



H 



I 



-^Bous h ■> f^uu + ^ hjjy . 



We begin with Type I critical points; for them Uc < 0. It is convenient to change 
the sign 

u I— > —u. 

Consider the string equation 

(5.5) S j [xu + t -Dbous^ - -Dbous/] dx^O 

for arbitrary two solutions f — f{u,v), h — h{u,v) to 

fuu — ^ fvv 1 U hyy . 
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Here the solution h = h{u,v) is related with the Hamiltonian (5.4) of the system (5.3) 
(with u replaced by —u) by 

(5.6) d^h = h. 

At the critical point {uc, Vc) one must have 

/° = /° = 

J UV V "'Cj vv 

J V 



J vv 

uvv 



+ fvv 



see (3.11 ) and (3.16) above). Let us expand the string equation near the critical point 



in the Taylor series with respect to the Riemann invariants coordinates 

2 

^ 3 

It is also convenient to perform a linear change of independent variables {x, t) 

(x+, xJ) 

X±=X+ ± v^/l°) t 



(cf. (3.17)). After the rescaling 
(5.7) 



of the shifted variables 

x± = x± — x\, f± = r± — r\ = V ±. a/u^ u 
at the limit A; — one arrives at the following system 



X. 



Z\/Uc 



(5.8) 
(5.9) 



with some coefficients Aq, Bq, A2, (see below). In the above equations we denote 

d 



dx 



± 



From the first equation one has 

I- 



4m, 
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Substitution in (5.9) yields 

(5.10) x_ = Aox+f^ + Bofi + e^A2 [{d^f^f + 2f_ 9^ f_] + e^A^ dtf^ + O (k^^^) 
where 



An 



Bn 



2/^2 27^/^0 8^3/2 



' A . . J VVV V C J WW ' J uv 



A') = Q ^/v^Bf) 



The substitution 



with the choice 



and 



72 

A4 = —UcBq. 

5 



x_ = a X 
x+ = 7T 
f_ =/3f/ 



a = -2''/'3'/'u'J''b'J' 



7 = - 



1/7 



29/735/7^1/7^^-153/7 



reduces the equation (5.10) to the standard form 



X = ?7T- 



lf/3 + l(^'V2f/f/")+^t/H+0(e^/^). 



The equation (5.8) takes the following form 



Returning to the original variables yields the following final form of the Universality 
Conjecture for the Type I critical behaviour for generic solutions to any equation to 
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the Boussinesq hierarchy 

r+ = rl+ ^ [x-x,+ {hi + it - tc)] 



! .6/7^, ' .4/7^, + ) 



(5.11^ 

r_=r_+/3e/ t/| ^^^^^^ , ^ ] + ^ {e ^ ) 



Here the function U{X, T) smooth for all real values X, T is determined from the Pf 



2 



equation, as it was explained above in Section 4.6 For the particular choice 



h{u, v) = V 

one obtains the critical behaviour of solutions to the Boussinesq equation itself. The 
solutions are parameterized by functions / = f{u, v) satisfying 

fuu ^ fvv 

Let us now proceed to considering the Type II critical behaviour. In this case it is 
convenient to use the original choice of signs for Boussinesq equation. At the critical 
point Uc > 0. Given a pair of solutions / = f{u,v), h = h{u,v) to Tricomi equation 



(5.1), consider the string equation of the form (5.5). At the Type II critical point one 



has 

f = f = f = f = f = 

The Riemann invariants 

3 

are complex conjugate. For the unperturbed system they are holomorphic resp. anti- 
holomorphic functions of the characteristic variables 

x± = X + {huTi \fu K)ti 



I.e. 



dx- dx- 



0. 



The Universality Conjecture for the Type II critical behaviour first formulated for the 
case of focusing NLS by T.Grava, C.Klein and the author [18j says that the holomor- 
phicity remains valid also for the leading term of asymptotic expansion of the perturbed 
solution. Moreover it was conjectured in [18j that this leading term is given in terms 



of the tritronquee solution to the Painleve-I equation (see above Section 4.6). Let us 
reproduce the arguments of [18] adapted for the case of Type II critical behaviour of 
solutions to Boussinesq equation (i.e., choosing h{u,v) = v). 
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As above we expand the string equations in Taylor series near the critical point. 
After the rescaling 

x±^ kx± 

of the shifted variables 

x± = {x — Xc) ± i \fUc it — tc) 
f± = r± - r±(uc,fc) =v±i 

one obtains at ^ 

(5.12) x+ = A(f2 -eVr+) + C(A;^/') 

^ = \{fuvv - ^ /D . ^ = 4^ v^- 

A complex conjugate equation holds valid for the variables a;„, r_. 
In the above equation 

Performing the same rescaling in the Boussinesq equation we find 

Ut = 

vt + u,u., = O {k'/^) . 

This is equivalent to 

So, at the leading order we may assume that f+ is an analytic function in The 
equation (5.12) then becomes 

(5.13) x+ = \ {fl - eV5+r+) + O {k^'^) . 
Choosing 



and doing a change of variables 



with 



k = e'/' 

x+ = a Z 
f+ = (3W 



a = (6A/i')^/^ 
A2 



we arrive at the needed equation 

Z = 6W^ -W" + (e^/s) 
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(see Remark 5.2 in [18] explaining the choice of the fifth root of the complex numbers). 
Returning to the original variables one obtains the Universality Conjecture for the 
Type II critical behaviour for equations of Boussinesq hierarchy in the form 



V — Vc — i u^J'^ {u — Uc 



(5.14) 



e 



where Wo{Z) must be the tritronquee solution to the Painleve-I equation. Indeed, it is 
easy to see that the line 

Z=-(x-Xc + (hl-iJu~ch°^(t-Q), xeR 
a 

belongs to the sector 

arg Z\ < —. 
5 



5.2. Toda hierarchy and its extension. The linear PDE (2.23) is of hyperbohc 
type for all real (m, v): 

(5.15) fuu — 6 fvv 

The D-operator transforming solutions of this equation into conserved quantities for 
(extended; see [7j) Toda hierarchy reads 



2 u 2 A \ ^ 2 

~ f ^ \_{fvv + 2 fuvv) "^x "I" 2 6 i^fuvv "I" 2 fvvv Ux ^x)] "I" ^ i 720 '^^^ 

+ Y^e" {ul + Uxx) + -^^fuvv [2 e" [u^ - 3ul^) + 3 u^^ vl + 24 vl^ ] 
+ [4 e" (2m^ Vx + 15 ul. v^x + 18 Uxx ^^xx) - 3 v^^ 

+ ^/.... [I6e2" (3wL-«') -4e" (7^^ t;^ + 8 t;^ - 12 1;^) - t^^] 
"2^""^^""=^''^ (I4e"n^ + I5t;^)-^A4. (iTe^^ + 27e« t;^ + 4 1;,^) 



52 



B. DUBROVIN 



The quasitriviality transformation is generated by the Hamiltonian 



3 



,2m 



2 2 

W ■ 

XX X 



^2u „A 



X X •^xx 



18 " uxx Vx -^e'^ul vl vl, - 6 u^x vl vl, + Avlvl, + 2 e''' ul vl, + 6 vl vl,] 

^3 

, (2 " + 14 « m,,. - 9 ^2 ^ ^4^, 



2880 A2 

c3 



C „ 4 C 

e u,v, — 



1152 A 



3840 



6 2 (U^ +4U3;3;j log 



f ^ — e 2 M 



where 



(cf. CT ). 

The main scheme explained above for the case of Boussinesq equation works es- 
sentially without changes for the Toda (and for the closely related NLS, see below) 
hierarchjj^ We will not enter into details here. 



5.3. Nonlinear Schrodinger equation. The D-operator for focusing NLS is given 

by 

J2 



f hf = Dnls/ = / - 



12 



fuuu ~l~ „ fuu 1 -\- 2 fy^u^UxVx ufuuu'^x 

lu 



(5.16) 



120 



Juuuu + 2^ fuuu ) U„ + 2U^^UxxVxx uf„U 



If 

csr^ JuUUu'UjxxV , JuUUvl^XxUx 

80 48m 



1 



3456^3 



(30/„„„ - 9ufuuuu + l^U^fbu + 4:U^f6u) ul 



^or> O ( ^fuUUV ~\~ ^UfuuUUV ~\~ 2lt f^Uv) Uj,Vx ~\~ „„„ (Q fuUUU ~l~ ^Ujf^u ~\~ 2lt ffju^ U, 

AiZw^ 288m 

1 M 1 

+ 2160 ^^^™ + lO^/s-) - ^ (18/5M + ^uhu) vl\+ 0(e6) 
For the defocusing case one has to replace u ^ ~u. In particular, taking 

f = \{uv^-u^) 



^x 



'^In the recent papers |35J the methods of [T7J [TB] were appHed to the description of the critical 
behaviour of certain particular solutions to the Toda hierarchy and some related systems. 
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one recovers the Hamiltonian of the focusing NLS equation 

1 

In this case the infinite series truncates. Taking g = —\v'^ + ^(logM — 1) one obtains 
the Hamiltonian of the Toda equation 

1 

hg = --v^ + u{\ogu - 1) - {ul + 2uvl) 



^' ' V 240^3 ^ 60m2 ^ 40m3 144m5 24^^ ^ 360^3 ' ^ ^ 

written in terms of the function = log u in the form 



e 



2, 



Remark 5.1. In [TB] a somewhat different scaling procedure was used for the descrip- 
tion of critical behaviour in the focusing NLS. 

5.4. Ablowitz - Ladik equation. Ablowitz - Ladik (AL) system [1] is written in the 
(complexified) form as follows 



(5.1^ 



ian = - ^{l- a„6n)(an-i + Ctn+l) + an 



ihn= ~ anhn){hn~l + hn+l) - hn. 



The reductions 

(5.19) hn = Tan 

correspond to the focusing/ defocusing versions of the AL equation respectively: 

(5.20) i hn + ^(a„+i - 2a„ + a„_i) ± |a„|^ ^""^^ ^"""^ = 0. 

The latter equation is universally recognized as a discrete version of the NLS equation. 



The system (5.18) is Hamiltonian with respect to the Poisson bracket 

(5.21) {a„, km} =i{l - a„6„) 6n,m, {an, am} = {&n, b,n} = 

with the Hamiltonian 

(5.22) = X] ^{anbn-i + 6„a„-i) + log(l - a„6n). 

n 

Actually the Hamiltonian 

(5.23) C = Y, log(l - «n&n) 

n 

generates phase transformations 

{fl-nj C} = an 

{bn,C}= ibn- 
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Factoring phase transformations out can be achieved by introducing new variables 
(5.24) 



log(l - a„6„) 



log 



log 



hn-l 



Observe that the variables Vn-, Wn take real values for the focusing/defocusing reductions 
of AL system. The variable Wn takes negative values for the focusing case and positive 
for the defocusing case. In the new coordinates the dynamics is governed by the 
Hamiltonian 



(5.25) 



// = ^ V(l-e-"'")(l-e 



COSVn 



with the Poisson bracket 

(5.26) {Wn, Vm} = Sn,m-1 " Sn,m, {w„, Wm} = {Vn, f m} = 0. 

After the interpolation and time rescaling 

Wn = w{en,et), Vn = v{en,et) 
one obtains, in the limit e — > the 1st order quasilinear system 

Wt = dx [(e~"' — l) sinf] 



(5.27) 



Vt = [e cos v~\ 



We call (5.27) the long wave limit of the AL system, in analogy with the long wave 



limit of Toda lattice. 



The characteristic velocities for (5.27) are 



A-t = — e ^ [sinv ± \^ — I cosf] 



So the long wave limit (5.27) of the AL system has elliptic type for the focusing case and 



hyperbolic type for the defocusing case. In the defocusing case the Riemann invariants 



of the system (5.27) are given by the expressions 
(5.28) 



r± = V ±2 arctan y/e"^ — 1. 



For a later convenience we will replace the dependent variable w{x) by a new variable 
u{x) by putting 



(5.29) 



u 



Ida: _ 1 



w 
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(cf. (2.8[)). The Poisson brackets of the variables u, v take the standard form (2.14) 



1 — exp{ 



1 - e^^^ 



The Hamihonian of the interpolated system reads 

Hal = j hAL dx = 
(5.30) 



u} \ ( 1 — exp{- 



1 



fiAL = (l - e ") cost; + — e '^u^ 



(e" — 4) cosf + 4(e'' — 1) sint> t>^ 



u} ) cosv dx 



In the long wave limit equations (5.27) one has just to replace the letter w with the 
letter u. 

Let us now establish a connection of the theory of AL equation with Hamiltonian 
perturbations of the nonlinear wave equation with the potential 

(5.31) P{u) = L«2(e-"). 

Recall that the dilogarithm is obtained by analytic continuation of the function 



k>l 



Li2{x) = ^ ^ for 1x1 < 1. 



The associated nonlinear wave equation reads 



(5.32) 



Utt 



- 1 



The linear PDE (2.23) for the conserved quantities of (5.32) reads 
(5.33) U = (e" - 1) fun. 



Lemma 5.2. The function f{u,v) is a conserved quantity for the long wave limit (5.27) 
of AL system iff it satisfies (5.33). 



Proof is given by a straightforward calculation. 



□ 



We can now apply the above arguments to constructing the D-operator for the AL 
hierarchy. 



Theorem 5.3. For any solution f = f{u, v) to (5.33) there exists a unique Hamil- 
tonian 



(5.34) 



Hf = Hf + Y,e'Hf= I DaU 

k>l 

Hf=jhfdx, hf=f 



dx 



commuting with the AL Hamiltonian (5.30).' 



6H „ SHal , SH ^ SHal 



6u{x) 6v{x) 6v{x) Su{x) 



dx = 0. 
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Explicitly 



1 



1 2(l-e«)/„„, + (e"-2)/, 



24 



vv 2 



, u vvv I fuvv fvv 2 
H ^iz'i^x H 

6 e«-l 12 e"-l 



The basis of linearly independent solutions to the linear PDE (5.33) can be obtained 
in the form 



(5.35) 



f^ = cosnve ""F(— n, — n; 1 — 2n; e") 
/~ = sinriw e~"'"F(— n, —n; 1 — 2n; e") 



n G Z, n 0. 



Here F(a,/3;7;x) is the Gauss hypergeometric function. For positive n the hyper- 
geometric factor is a polynomial in e" of the degree n; these conserved quantities 
correspond to the first integrals of AL hierarchy constructed in The second part of 
the solutions to (5.33), the one for n < 0, corresponds to the new integrals of the AL 



hierarchy found recently in 
For example, choosing 

one obtains 

The choice 

yields the e expansion 



/ = /+ = (l-e-") cos^; 



hf = hAL- 



f = fr 



'1 — e ") sinf 



;i -e sini; + — e X- 



(g« _ 4^ gin vUx — 4(e" — 1) cos v 
e" - 1 



of the Hamiltonian 



of the discrete modified KdV equation 

1 



n — 2 (-'- ~ 0"nbn){0'n+\ — O'n-l] 



hn = - {I- a„6„)(6„+i - hn-i] 



Observe that the Hamiltonian of the perturbed wave equation (5.32) reads 

ul - (e« - 1) {ul + 2vl 



(5.36) 



pert 



-t;2 + Li2(e-") + 
2 ^ ^ 24 



(e" - 1)^ 



dx + O (e^) . 



This perturbation is not equivalent to the one given in (4.19) coming from the gener- 



alized FPU system with the potential (5.31) 
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Remark 5.4. The recurrent procedure of constructing the conserved quantities of the 
long wave limit of the Ablowitz - Ladik hierarchy can be described in terms of the 
following solution to WDVV 

(5.37) v) = ^uv^- Li3 (e"") 



(see Remark 2.1^ above). 
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